X,Y,Z- Waves: Extended Structures in Nonlinear Lattices 
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Motivated by recent experimental and theoretical results on optical X-waves, we propose a new 
type of waveforms in 2D and 3D discrete media - multi-legged extended nonlinear structures (ENS), 
built as arrays of lattice solitons [tiles or stones, in the 2D and 3D cases, respectively). First, we 
study the stability of the tiles and stones analytically, and then extend them numerically to complete 
ENS forms for both 2D and 3D lattices. The predicted patterns are relevant to a variety of physical 
settings, such as Bose-Einstein condensates in deep optical lattices, lattices built of microresonators, 
photorefractive crystals with optically induced lattices (in the 2D case) and others. 
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Introduction and setup. In the past few years, stud- 
ies of dynamics in discrete media, as well as in their 
continuum periodically modulated counterparts, have en- 
joyed very rapid progress This is fueled by a need to 
tackle diverse physical contexts such as nonlinear optics 
of waveguide arrays and photorefractive crystals j2i]j the 
physics of Bose-Einstein condensates (BECs) in optical 
lattices the denaturation of the DNA double-strand 
, and so on. These efforts have resulted in the creation 
of a variety of novel nonlinear structures, such as discrete 
dipole 0, quadrupole necklace f?! and other multi- 
ulse/multi-pole localized patterns iS], discrete vortices 
ring solitons and others ^J. Such structures 
have a potential to be used as carriers and conduits for 
data transmission and processing, in the context of all- 
optical and quantum- information schemes. 

The present work is motivated by the above-mentioned 
achievements, as well as by recent developments in stud- 
ies of X-waves, which are quasi-linear extended structures 
generated in three-dimensions (3D) |12|. two-dimensions 
(2D) 13], and even quasi-discrete [1J| media. These 
waves arise when the linear second-order differential 
operator in the relevant wave equation is genuinely, 
or effectively sign-indefinite (roughly speaking, a 

D'Alembertian instead of the Laplacian). A similar shape 
is featured by cross configurations formed by two in- 
tersecting domain walls in a system of two coupled 2D 
Gross-Pitaevskii equations with a difference that 

these structure are essentially nonlinear ones. 

In the present work, we are dealing with the sign- 
definite discrete Laplacian (nevertheless, one can easily 
induce the sign-indefiniteness in the discrete wave equa- 
tion, applying the well-known staggering transformation 
Pi along one coordinate only). Instead of emulating the 
mathematical structure of the continuous wave equations 
supporting X-waves in the lattice medium, our aim is to 
construct stable extended nonlinear structures (ENSs) 
in discrete media. These wave structures are partly de- 
localized, by their nature, ranging from a few sites to 
infinitely extended ones. An essential difference between 
the quasi-linear X-waves in continua and the lattice enti- 
ties presented here is that the latter are strongly nonlin- 
ear structures, being built as arrays of individual lattice 



solitons ("tiles"); in that sense, t hey somewhat resem- 
ble weakly localized hypersolitons |l7| and supervortices 
[Tsj l. constructed as finite ensembles of individual soli- 
tons. In fact, we initialize our analysis (see below) from a 
purely nonlinear (anti- continuum, AC) limit, without lin- 
ear intersite coupling on the lattice. The eventual shape 
of the extended lattice waves features a great variety, 
including X-waves, Y- waves, Z- waves and other forms, 
provided that the underlying core ( "tile" ) is dynamically 
stable. Each species is stable in a finite interval of values 
of the intersite coupling constant. Some of those pat- 
terns may definitely be of interest (for instance, as star- 
shaped waveguiding configurations) to applications in 
both nonlinear-optical and matter-wave contexts, where 
they can potentially be created. 

The analysis is performed in terms of the discrete non- 
linear Schrodingcr (DNLS) equation in 2D or 3D, 
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where Un is a complex amplitude of the electromagnetic 
wave in an optical waveguide array (in the 2D case) 0, 
[T^, or the BEG wave function at nodes of a deep (2D 
or 3D) optical lattice [l^, n being the vectorial lattice 
index, and A2 the standard discrete Laplacian. Further, 
e is the constant of the intersite coupling, and the overdot 
stands for the derivative with respect to the evolution 
variable; the latter is z in optical arrays, or t in the BEG 
model (or in a crystal built of microresonators l2ll|). 

The presentation of our results on ENS solutions of Eq. 
proceeds as follows. First, in the AG limit (e = 0), 
we specify main types of the "tiles" , of which the ENSs 
are to be composed. Then, we numerically extend these 
compositions into full ENS solutions for e > 0. Finally, 
we summarize our findings and discuss further problems. 

Analytical results: tiles, stones, and their stability. We 
look for standing- wave solutions to Eq. JQ), in the form 
of Un = exp(iAt)(/)„, with 0n satisfying the equation. 



/(0„, e) = A0„ - eA20„ - |0„p0„ ^ 0. 
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Perturbation of Eq. around the solutions of Eq. 
leads to the linearization operator. 
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the overbar standing for complex conjugate. By means of 
a rescaling, we fix A = 1, keeping e as a free parameter. 
In the AC Hmit of e = 0, solutions to Eq. (O are 
= re*^" , where real amplitude r is or -s/A, and 9^ are 
arbitrary constants. To continue such a solution to e ^ 
0, the Lyapunov-Schmidt condition needs to be satisfied 
|2^ . viz., the projection of eigenvectors of the kernel of 
Hn^ (i.e., zero modes of the operator) onto the system 
of stationary equations should be null. This condition 
gives rise to a solvability condition at each "AC-fiUed" 
site, i.e., one with r in the AC limit. 
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(the factor of —2i is introduced for convenience). Impor- 
tantly, the eigenvalues 7 of the Jacobian, Mij = dgt/dOj, 
are intimately related to the leading-order approxima- 
tion for eigenvalues (A) of the linearization around a 
stationary solution, A = ±-\/27. We use this rela- 
tion, alongside a perturbative expansion in the solution, 
4>n = 4>ri' +e<^n ^ + • ■ ., to dcrivc leading-order bifurcation 
conditions for any given configuration, and subsequently 
to find the corresponding linear-stability eigenvalues. 

2D configurations can be categorized by the number of 
their "legs" (i.e., the number of quasi-lD lines of excited 
sites it contains). Simplest is the one-leg structure as- 
sembled of two- or three-site tiles, for which ID stability 
results are applicable here, the resulting Jacobian is 
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from which conclusions for the stability of ENSs may be 
drawn. In particular, a fundamental stability condition 
extensively used in our analysis is that the phase shift be- 
tween fields at adjacent sites must be tt |23jI23- -F'or the 
two-site tile satisfying this condition, the stability eigen- 
values are A ~ ^/2rj ~ ±2-/ei, while for its three- 
site counterpart, (+1, — 1,+1), they are A — ±\/2e« and 
A = ±-\/6e«. The full one-leg configuration assembled of 
these tiles is a chain, (. . . , —1, -1-1, —1, -1-1, —1, . . .), here- 
after termed configuration 1. 

Proceeding to two-tile configurations allows us to ex- 
amine more complex (in particular, X- and Y-like) 2D 
structures. X-configurations can be built of square- or 
cross-shaped tiles, which, in the AC limit, are 




and 




(5) 



However, the latter one gives rise to an unstable eigen- 
value pair, A — ±2e (to leading order), therefore it is 
not considered further. The former tile has a pair of 
A = ±2\/2ei, and a double one, A ±2ei, hence stable 
X-configurations (to be called 2a) can be constructed at 
finite e as arrays of such tiles with alternating signs. 
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FIG. 1: (Color online) A family of single-leg patterns in the 
2D DNLS equation. The top left and right panels show, re- 
spectively, the norm and the largest instability growth rate vs. 
e. The contour plots in the middle panels show the stationary 
solutions for e = 0.05 (stable) and e = 0.1 (unstable). The 
spectral planes (Ar,Ai) of numerically computed eigenvalues 
for these solutions are displayed in the bottom panels. 




FIG. 2: Same as in Fig. 1 but for configurations 2a, 2b and 2c. 
Bottom panels present the stability eigenvalues for individual 
tiles of which these patterns are built, numerically found ones 
(solid lines) versus the analytical prediction (dashed lines). 
Examples of species 2a and 2b are displayed in two left and 
two middle panels for e = 0.15 and 0.25, and an example of 
2c in the right panels pertain to e = 0.05 and e = 0.1. 
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FIG. 3: Same as Fig. 2, but for configurations 3a and 3b. 
Examples of both species are shown for e = 0.05 and 0.1. 



Other two-leg configurations are Y-shaped and 
X" ones, composed, respectively, of the tiles: 



and 





'skew- 
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Both these tiles are stable, with respective eigenval- 
ues A — ±^/2ei and A = iy/6ei for the former, and 
A = ±0.8747ii, A = ±lM3y/ei, A = ±2.882^^^ and 
A = ±2.690-yei for the latter. Stable sign-alternate ENSs 
assembled of them are termed 2b and 2c, respectively. 

Three-leg ENSs can be built of other (generally, denser 
filled) stable tiles, such as 
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The former tile possesses three pairs of double eigenval- 
ues, ±^/2e^, zb-s/Oei, and ±^/8ei, and two pairs of single 
ones, ±.2y/li and ±Vl2ei. The latter has a double eigen- 
value ±-\/2ei and a single one ±-\/8e«- Three- leg config- 
urations constructed of these tiles with alternating signs 
are called 3a and 3b, see below. In fact, the latter one 
will be a Z-shaped array. 

Finally, we consider two configurations as a proof-of- 
principle of the extension of the ENS concept to the 3D 
space. In particular, augmenting the first tile of by 
two —1 sites, adjacent to middle 1 along the third direc- 
tion, we fabricate a "stone" . The stone is stable, with sin- 
gle eigenvalues A = ±1.248-^61, A = ±2^/ei, A = ±^/6€i, 
A = ±2.763^/ei and A = ±3.848-v/e«, and double and 









I 


*■■> 


1 ' h 


1 


1 


1 


















1 




* 


t 






il 






il 


1 

1 




















1 






•I 


1 




% : '= 




-1 


1 


-1 







1 1 1 lis 11 1 1 I i 11 
'l I 



41 1 1 Hi 11 41 415 1 Hi 11 
I, I, 





0.01 0.02 E 0.03 0.04 0.05 



FIG. 4: Same as in above figures, but for 3D configurations 4a 
and 4b. Their examples are plotted, in two left and two right 
panels, respectively, for e = 0.03 and e = 0.1. In configuration 
4b, the imaginary parts are shown by different colors, green 
for phase 7r/2 and yellow for 3n/2. 



triple ones, A = ±^/8ei and A — ±V2ei, respectively. 
These stones will be used to build a stable 3D pattern 
(again, with sign alternations), called 4a in our nomen- 
clature. 

The simplest stable stone has a "diamond" -like shape, 
i.e., it consists of six AC-filled sites surrounding, as near- 
est neighbors, an empty one. Although some configura- 
tions of this type are unstable ^ , we have found a stable 
one, carrying a phase distribution that corresponds to a 
quadrupole in the plane, with phases 7r/2 and 37r/2 lent 
to the two out-of-plane sites. This stone has triple and 
single eigenvalues, A = ±2ei and A = ±4ei, respectively. 
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and a higher-order one, C(e^). Augmenting stones of this 
type by alternating phase pulses, we will build a stable 
3D structure labeled 4b. 

Numerical results. In Figs. 1 1141 we present, in a unified 
format, results of the numerical continuation of various 
ENS configurations from the AC limit. Figure shows 
configuration 1 (a straight chain of sign-alternating tiles) . 
The norm of the configuration, N = (j)^, and its max- 
imum instability growth rate are shown, as a function of 
e, in the top left and right panel, respectively. In the 
21x21 lattice used in our 2D computation, the configura- 
tion becomes unstable for e > 0.074. The middle panels 
show typical examples of stable and unstable configu- 
rations, and the bottom panels show the spectral plane 
(Ar,Ai) of their numerically found (in)stability eigenval- 
ues, A = + the instability being ushered by Ar 7^ 0. 

Figure|2]presents three different varieties of the two-leg 
configurations: 2a (X-waves, left panels), 2b (Y-waves, 
middle panels) and 2c (skewed X-waves, right panels). 
In addition, numerically computed stability eigenvalues 
of the corresponding tiles, i.e., the first one in Eq. Q 
and both in Eq. 10, are shown by solid blue lines for 
comparison with the analytical predictions of the previ- 
ous section (dashed red lines), indicating good agreement 
between them. Configurations 2a, 2b, and 2c become un- 
stable at e > 0.202, e > 0.206 and e > 0.068 respectively. 

Similarly, Fig. Oshows two three-leg configurations, 3a 
and 3b (the latter one may be naturally called a Z-wave), 
and their stability characteristics. These configurations 



become unstable at e > 0.053 and e > 0.075, respectively. 

Finally, figure 0] shows the continuation to e ^ of 
3D configurations 4a and 4b (defined above), which are 
stable, respectively, at e < 0.043 and e < 0.075. Note 
that the corresponding numerically computed eigenval- 
ues (for the respective stones) again agree well with the 
analytical predictions. 

Conclusion. We have presented a systematic approach 
towards constructing a variety of extended states on 2D 
and 3D nonlinear lattices. They are composed of building 
blocks, namely, tiles (2D) or stones (3D) with alternating 
signs, that are originally defined in the anti-continuum 
limit (the sign alternation is necessary for stability of 
the patterns). Stability intervals for the structures were 
revealed by numerical continuation of the corresponding 
solution families to finite values of the coupling constant. 
Examination of nonlinear evolution of unstable patterns 
is another relevant problem, to be considered elsewhere. 

Another interesting topic is a discrete model such as 
Eq. (^), but with a negative intersite coupling in one 
spatial direction (say, x); as mentioned above, this model 
can be generated from Eq. by means of the staggering 
transformation applied in this direction, i.e., u„^.„^.„^ = 
[—l)'^'''Vn^.ny,n^- In the 2D version of such a model, one 
can immediately find two exact linear single-leg solutions 
oriented along the diagonals, viz., Un^,n = exp(— jria; =p 
Uyl). Their superposition can be used to construct exact 
linear X-waves. A potential extension of the latter states 
into the nonlinear model would also be worth considering. 
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